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, We interpret hadronic particle abundances produced in S-Au/W/Pb 200 A GeV reactions in terms of 

the final state hadronic phase space model and determine by a data fit of the chemical hadron freeze-out 
parameters. Allowing for the flavor abundance non-equilibrium a highly significant fit to experimental 
'sj" ' particle abundance data emerges, which supports possibility of strangeness distillation. We find 

^ , under different strategies stable values for freeze-out temperature T{ = 143 ± 3 MeV, baryochemical 

potential /ie — 173 ± 6 MeV, ratio of strangeness (7s) and light quark (7q) phase space occupancies 
7s/7q = 0.60 ± 0.02, and 7q = 1.22 ± 0.05 without accounting for collective expansion (radial flow). 
' When introducing flow effects which allow a consistent description of the transverse mass particle 

O '. spectra, yielding |wc| = 0.49 ± 0.01c, we find 7s/7q = 0.69 ±0.03, 7q = 1.41 ±0.08. The strange quark 

00 ■ fugacity is fitted at As — 1.00±0.02 suggesting chemical freeze-out directly from the deconfined phase. 

1 Introduction 

In relativistic nuclear collisions, for a short period of time, physical conditions are recreated similar to 
. those present in the early Universe about 40 iis after the Big-Bang. The primary objective of the extensive 
experimental program at BNL (Brookhaven National Laboratory) and CERN (European Organization 
for Nuclear Research) is the understanding of 'vacuum melting', the formation and properties of the 
deconfined quark-gluon plasma (QGP) phase of strongly interacting (hadronic) matter 
5^ I Our present investigation addresses the completed experimental analysis of the CERN 200 A GeV 

Sulphur beam reactions with laboratory stationary 'heavy' targets, such as Gold, Tungsten or Lead nuclei 
This highest presently available energy content materializes in form of hadronic particle multiplicity and 
corresponds to nuclear collisions occurring in the center of momentum frame Eqm = V^/B = 8.8 GeV = 
9.4mNC^ for each participating nucleon. Several prior attempts have been made to interpret the hadronic 
particle abundances and spectra within statistical particle production interpretation ^, |^, |^, |7[ 
and the present report takes, as we believe it, to a conclusion this extensive body of research work. 
The result of our present analysis are the parameters and properties of the hadronic system at chemical 
freeze-out. The new theoretical element which makes this attempt successful and which we introduce 
here is the chemical non-equilibrium abundance of light quark flavors. 

This analysis of the experimental data is based on the assumption that a relatively small and dense 
volume of highly excited hadronic matter, the 'fireball', is formed in the reaction. When we speak of 
particle abundance freeze-out, we refer to the stage in evolution of the fireball at which density has 
dropped to the level that in subsequent collisions particle abundances remain unchanged. We side-step 
the need to understand largely unknown collective flows originating in both, the memory of the initial 
'longitudinal' collision momentum, and the explosive 'transverse' disintegration driven by the internal 
pressure of the hadronic matter. Thus, all particles originate from an expanding fireball surface. In order 
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to minimize the dependence on the source shape, geometry and dynamics, we either consider ratios of 
particles obtained at central rapidity, or ratios of 47r-abundances. Suggestions were made that detailed 
production dynamics needs to be considered [§, allowing that some particles are emitted well ahead 
of others, and of the final state hadronization. In our approach we see evidence for at best, two stage 
process, with some distillation of flavor content, and specifically strangeness leading to s-quark enriched 
final stage. However, there is no compelling evidence for a real multi-stage, continuous particle production 
process in our evaluation of the data. 



2 Chemical non-equilibrium statistical model 

Along with the other statistical analysis of the experimental data mentioned above, we employ the 
local thermal equilibrium method, and use a local freeze-out temperature Tf. One can argue that the 
accessibility of many degrees of freedom, as expressed by high specific entropy content, validates use 
of the thermal equilibrium. The usefulness of the local thermal equilibrium is, in part, based on the 



experimental fact that, at given transverse mass m± = ^Jm? ~ 1.5-2.5 GeV, mj_-spectra of very 
different hadronic particles have the same inverse slope Tq |1C]. To a very adequate precision, for these 



high m^, this inverse slope Tq ~ 235 MeV arises from Tf according to the Doppler-blue-shift relation 
using a local collective flow velocity. Impact of flow effects is discussed in detail in section ^. 

Regarding chemical equilibration, we will recognize here two different types, relative and absolute 



[12 1, illustrated in the following example. Consider within a fixed volume a hot, thermally equilibrated 
gas of nucleons N = p,n, pions 7r^,7r^, and A-resonances: 

• Relative chemical equilibrium 

The relative abundance of the u and d carrying quarks is easily established through flavor exchange 
reactions, such as p + 7r~ <-> n + n^; at the valence quark level there is no creation or annihilation process 
that need to occur in each reaction to equilibrate, between the different hadronic particles, the quark 
flavors. We speak of relative chemical equilibration occurring through quark exchange reactions described 
in terms of fugacities Xi, i = u, d, s. Prom quark fugacities, the particle fugacities are reconstituted, and 
in particular we recall that baryochemical potential arises as = 3TlnAq, where we assumed that the 
difference between u and d quarks is not significantly affecting the particle ratios studied here |p. 

• Absolute chemical equilibrium 

Very much different is the process involving the equilibration between the number densities of mesons and 
baryons. Baryon-antibaryon formation processes are of the type N + N <^ p -\- oj and require that aside 
of the reorganization of the quark content, also the number of valance quarks changes, here between 6 
and 4. The valance quark numbers has to be controlled by new parameters, the phase space occupancies 
7j, i = u, d, s. Only when 7j is present can the freeze-out abundance of hadron families be differentiated. 

The difference between Aj and is that, e.g., for strange and anti-strange quarks the same factor 7s 
applies, while the antiparticle fugacity is inverse of the particle fugacity. The proper statistical physics 
foundation of ji is obtained considering the maximum entropy principle for evolution of physical systems. 
In such a study it has been determined that while the limit 7^ ^ 1 maximizes the specific chemical entropy, 
this maximum is extremely shallow, indicating that a system with dynamically evolving volume will in 
general find more effective paths to increase entropy, than offered by the establishment of the absolute 
chemical equilibrium |11|. 



Generally, the microscopic processes which lead to absolute chemical equilibrium are slower than 
those leading to establishment of thermal or relative chemical equilibrium, considering that the inelastic 
reactions that change the particle abundance have usually much smaller cross section than particle number 
preserving exchange reaction, or elastic collisions. While thermal, and relative chemical equilibrium will 
in general occur within the life span of dense matter formed in nuclear collisions, approach to absolute 
chemical equilibrium provides interesting chronometric and structure information about the dynamics 
of the collision process |jl^, ^]. In consequence, the number density of different hadrons produced in 
the collision has to be characterized by individual fugacities. If each hadronic species were to require its 
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own fugacity, there would be little opportunity to relate reliably the observed particle yields to specific 
physical processes. 

The interesting insight of this paper, concluded from the success of the description of the experimental 
data, is that chemical non-equilibrium is important for the understanding of the experimental results, 
and that there is need for just four chemical abundance parameters: two standard fugacities Aq and Ag 
(relative chemical equilibrium), and two chemical non-equilibrium occupancy parameters, measuring the 
abundance of valance quarks, 7s, a parameter describing the approach of strange quarks to phase space 
abundance equilibrium ||2|, 13], and 7q, introduced here for the first time into such data analysis. 



characterizing the phase space abundance of light quarks. 7q is required as a parameter in the picture of 
direct hadronization of a deconfined QGP region into final state hadrons, which are free streaming and do 
not form an intermediate hadronic phase. When developing a dynamical description of QGP evaporation, 
one immediately establishes that the entropy excess we reported previously [14|, requires a hadronic phase 



space parameter capable to control the excess of hadronic particle abundances consistently. Moreover, 
considering that in the fragmentation of gluons additional light quark pairs arise, and can over populate 
the phase space abundance of light quarks, we are naturally lead to implement the light quark occupancy 
parameter 7q. Note that such direct emission from QGP scenario, while convenient to understand the 
meaning of our data analysis in microscopic terms, is by far not a unique framework for the interpretation 
of the results we obtain. Which microscopic model is appropriate will be understood when other systems 
can be subject to a similar analysis as the one presented here. Comparable particle abundance data is 
presently assembled both for Au-Au collisions at 11 A GeV at AGS-BNL, and for Pb-Pb collisions at 
158 A GeV at CERN, with intermediate energies available at CERN in near future. 

The evaluation of the final particle yields follows the pattern established in earlier work (see, e.g., 
III, ^). The relative number of primary particles freezing out from a source is obtained noting that 
the fugacity and phase space occupancy of a composite hadronic particle is expressed by its constituents 
and that the probability to find all j-components contained within the i-th emitted particle is: 

N,^Y[j,X,e-^^/^, A, = nA,, l^ = Il^J■ (1) 

jet j& j& 

The unstable hadronic resonances are allowed to disintegrate and feed the stable hadron spectra. Central 



rapidity region y ~ 0, where Ei = ym? + p"^ = \Jmf +p^coshy , or the full phase space coverage as 
required by the kinematic range of the experiments are considered. Since > 1 GeV hadron spectra are 
not significantly deformed by decays (see Figs. 3 and 4 in Ref. [P), one may in this limit first evaluate the 
partial multiplicities for different hadrons, and than allow these to decay, which speeds considerably the 
calculations. We have verified that this approach, which is exact for Att multiplicities, gives results which 
are much more precise than the experimental uncertainties for the central y ~ high p± data, and thus we 
have used this approach. Once the parameters Tf, Aq, Ag, 7q, 7s are determined from the particle yields 
available, we can reconstitute the entire hadronic particle phase space and obtain the physical properties 
of the system, such as, e.g., energy and entropy per baryon, strangeness content. Even though we are 
describing a free streaming gas of emitted particles, we can proceed as if we were evaluating partition 
function of system with the phase space distributions described by the statistical parameters, given that 
in a gedanken experiment we have just in an earlier moment still a cohesive, interacting system. We 
have implemented all relevant hadronic states and resonances in this approach and have also included 
quantum statistical corrections, allowing for first Bose and Fermi distribution corrections in the hadron 
abundances and in the phase space content. These corrections infiuence favorably the quality of the data 
fits. 

Most of our analysis will be carried out within this thermo-chemical framework. However, in section |5| 
we shall enlarge the discussion to include the analysis of m±_ particle spectra, introducing the collective 



velocity Vc. The different schemes that are possible to implement fiow were studied |15|. We adopt here 
a radial expansion model and consider the causally disconnected domains of the dense matter fireball 
to be synchronized at the instance of collision — in other words the time of freeze-out is for all volume 
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elements a common constant time in the CM frame. The freeze-out occurs at the surface of the fireball 
simultaneously in the CM frame, but not necessarily within a short instant of CM-time. 

Within this approach the spectra and thus also multiplicities of particles emitted are obtained replac- 
ing the Boltzmann factor in Eq. (|l|) by []l5|| : 

e-^^/^ - ^ / dn^j.il + V, ■ p,/Ej)e-"-^^'+-^-P^/'^^^ , 7. = , (2) 

a result which can be intuitively obtained by a Lorentz transformation between an observer on the surface 
of the fireball, and one at rest in laboratory frame. In certain details the results we obtain confirm the 
applicability of this simple approach. 

As can be seen in Eq. (^) we need to carry out an additional two dimensional (half-sphere) surface 
integral with the coordinate system fixed by an arbitrary, but fixed (collision) axis which defines the 
transverse particle momentum. Just a one dimensional numerical integration over one of the surface 
angles needs to be carried out. To obtain the particle spectra as function of 'm± we need also to integrate 
over rapidity y . This rapidity integration can be approximated for a narrow rapidity interval using the 
error function. 

3 Analysis of experimental results 

We now discuss the results of our data analysis. In table |l], we present 10 different least square fits to 
the data obtained using MINUIT96.03 program from the CERN Fortran library. The two fits (A, A') are 
carried out allowing just two parameters, Tf and Aq, leaving the remaining three chemical parameters 
at their implicit value (= 1). The difference between the two fits is that A (and B,C,D) exclude data 
comprising the J7-yields for reasons which we will discuss momentarily. We consider here 18 data points 
listed in table § (of which three comprise the O's) which we fit with these two parameters, thus in fit A' 
we have 16 and in fit A 13 degrees of freedom (dof), while in the all-parameter fits D and D' we have 
13 and respectively 10 dof. As seen by the large x^/dof shown in last column of table |l|, the fit A is not 
acceptable. Moreover, fitting Ag (fits B,B'), and thus in effect testing the so called hadronic gas model 

is actually not improving the validity of the description, even though we have not enforced here the 
strangeness conservation among the emitted hadrons. Setting {s — s) = introduces a constraint between 
the 3 parameters which is difficult to satisfy. Implementing such strangeness conservation (not shown in 
table H) we find that the fit to all data has x^/dof ~ 24 , with 16 dof. 

In short, the equilibrium hadron gas model of particle production fails today to describe the precise 
experimental particle abundances, our finding contradicts strongly earlier expectations [^]. The major 
discrepancies, as shown in table |2| column B, include S/A (10 standard deviations (sd) and 7sd), H/A 
(6sd), i^s/A (5sd). Two values of same particle ratios appear here since if a particle ratio is available 
at fixed p_L and fixed m±, both are fitted and the total included in analysis, as these experimental 
results correspond to different data samples. We did not include in our analysis the same data sample as 
seen in Q, but the number of degrees of freedom remains similar. We have introduced precision central 
rapidity strange baryon and antibaryon results, and omitted several results (see below) that we could not 
consider as having a straight-forward thermal interpretation. We do not see our data selection as biased, 
as in the end we are retaining a quite diverse sample typically with smallest errors. We did not consider: 

• the NA44 deuteron yields since it is difficult to distinguish the direct, thermally produced fraction from 
final state interaction formation (and depletion); 

• rj/ir^ of WA80 as the fiavor evolution of the r] has not been fully understood within the thermal model; 

• (f)/ (p+io) of NA34 (Helios 3) as it is obtained in forward rapidity domain and thus subject to considerable 
fiow influence not studied in detail here; 

• the early NA44 pion to nucleon ratios were used in but these are obtained from non-overlapping 
rapidity regions, and data was not corrected for the weak hyperon feed down decay. However, we 
considered that it is very important to include a similar data sample in the fit, as the pion to nucleon ratio 
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Table 1: Statistical parameters obtained from fits of S-Au/W/Pb data. In fits A to D, particle abundance 
ratios comprising Q are not fitted. In fits A' to D' all experimental data in table |2| were used; in fits 
Dg and Dg strangeness conservation in the fitted particle yields was enforced. Asterisk (*) means a fixed 
(input) value, not a parameter of the fit. 



Fits 


Tf [MeV] 




■> 




Ag 






IS 




/q 


/dof 


A' 


145 ± 3 


1.52 


± 


0.02 


1* 






1* 




1* 


17 


A 


145 ± 3 


1.52 


± 


0.02 


1* 






1* 




1* 


21 


B' 


144 ± 2 


1.52 


± 


0.02 


0.97 ± 


0.02 




1* 




1* 


18 


B 


144 ± 2 


1.53 


± 


0.02 


0.97 ± 


0.02 




1* 




1* 


22 


C 


147 ± 2 


1.48 


± 


0.02 


1.01 ± 


0.02 


0.62 


± 


0.02 


1* 


2.4 


C 


147 ± 2 


1.49 


± 


0.02 


1.01 ± 


0.02 


0.62 




0.02 


1* 


2.7 


D' 


144 ± 3 


1.49 


± 


0.02 


1.00 ± 


0.02 


0.73 




0.02 


1.22 ± 0.06 


0.90 


D 


143 ± 3 


1.50 


± 


0.02 


1.00 ± 


0.02 


0.73 




0.02 


1.22 ± 0.06 


0.65 


D'. 


153 ± 3 


1.42 


± 


0.02 


1.10 ± 


0.02 


0.70 




0.02 


1.26 ± 0.06 


3.04 


D. 


153 ± 3 


1.42 


± 


0.02 


1.10 ± 


0.02 


0.70 




0.02 


1.26 ± 0.06 


3.47 



is strongly temperature dependent (see below Fig.|^. However, the results from NA35 are available only 
for S-Ag reactions, a somewhat lighter collisions system. We have looked carefully at the systematics of 
these and other results offered by the NA35/49 collaboration and have determined that the ratio /{p—p) 
{h^ are all positively /negatively charged hadrons) shows little variability between the different systems 
studied. In S-S the result is 4.6 ± 0.4, in S-Ag (see table |) it is 4.3 ± 0.3 and in Pb-Pb it is 4.4 ± 0.5, 
see Ip for data review. We thus adopted in our fit the central value observed in S-Ag collisions. The 
smallest error we chose is a conservative assumption here, as it in principle pushes the overall error of 
the fit up. In passing we also note that another related variable we consider here in difference to earlier 
work is the ratio {h~^ — h~)/{h~^ + h~). We fit now the An data of EMU05, rather than the central 
pseudo-rapidity domain value. 

In the fits (C, C), we introduce 7s, i.e., strangeness chemical non-equilibrium is allowed for Q. This 
single parameter reduces greatly x^/dof, which drops by a factor 10, providing a major improvement 
in comparison between theory and experiment. The resulting fit is fully compatible with the scenario 
B in [^. However, presence of additional and/or more precise experimental data makes this fit, which 
was 5 years ago possible, now statistically also not acceptable. For 18 — 4 = 14 dof (fit C including 
O's), x^/dof = 2.4 is expected to arise in less than 1% of experiments. Clearly some physics is missing 
in this description of the data. It has been argued that a more refined dynamical particle production 
model is needed to account for this discrepancy, refiecting on some complex dynamics involved in particle 
production (see, e.g., Q). The surprise is that all that is needed is the chemical non-equilibrium of the 
light quark abundance. 

We now introduce 7q. The full fit to all data (fit D') is statistically significant, with x^/dof = 0.9 . 
The profile of x^/dof as function of Tf shown in Fig.[l| (thin solid line) shows a well defined minimum 
at 144 lb 3 MeV. We note, in passing, that quite different considerations have lead others to propose 
a universal chemical freeze-out temperature at nearly exactly this value p^ . The thick solid line in 
Fig.|l] (fit D) arises when, in the fit, one excludes three data points involving completely strange O's. 
We came to try this, since when exploring the stability of the fit D' against suppression of some of the 
data, we noted that the fit D as expressed by x^/dof is better. As can be seen in Fig.|l], the fit D 
minimum (thick solid line at Tf = 143 MeV) is sharper and x^/dof drops to 0.65 . In the table |l|, we 
further see that the removal of these three data points has for all the other fits the effect of reducing the 
significance of the fit, as would be expected given the few degrees of freedom we have at our disposal. 
That an improvement occurs when the fit is statistically significant signals that there is probably an 
additional mechanism of 0, formation contributing to the yields. This indeed is what we expect, since 
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Table 2: Available particle ratios: experimental results for S-W/Pb/Au, references for the results in sec- 
ond column, momentum or transverse mass cuts (in GeV) in third column, followed by columns showing 
the different fits, corresponding to those in table |l[ Asterisk * means a predicted result (corresponding 
data is not fitted or not available). 



Ratios 


Rcfs. 


Cuts 


Data 




Fit A 


Fit B 


Fit C 


Fit D 


Fit D' 


Fit T)s 


Fit D 


S/A 


1 


1.2 <p+<i 


0.097 ± 





006 


0.162 


0.157 


0.105 


0.099 


0.100 


0.110 


0.111 


S/A 


1 


1.2 < Pi < 3 


0.23 ± 





02 


0.36 


0.38 


0.23 


0.22 


0.22 


0.18 


0.18 


A/A 


1 


1.2 < Pi < 3 


0.196 ± 





Oil 


0.194 


0.202 


0.205 


0.203 


0.203 


0.203 


0.203 




1 


1.2 < J3± < 3 


0.47 ± 





06 


0.43 


0.48 


0.44 


0.45 


0.44 


0.33 


0.33 




2 


p+ > 1.6 


0.57 ± 





41 


1.00* 


1.18* 


0.96* 


1.01* 


0.98 


0.55* 


0.55 


n+n 


2 


p+ > 1.6 


0.80 ± 





40 


0.27* 


0.27* 


0.17* 


0.16* 


0.16 


0.16* 


0.16 


k+iK- 


1 


p+ > 0.9 


1.67 ± 





15 


1.96 


2.06 


1.78 


1.82 


1.80 


1.43 


1.43 


K°/A 


3 


p+ > 1 


1.43 ± 





10 


1.51 


1.56 


1.64 


1.41 


1.41 


1.25 


1.25 


K^/A 


3 


p+ > 1 


6.45 ± 





61 


7.85 


7.79 


8.02 


6.96 


6.96 


6.18 


6.17 


K°IA 


\ 


m+ > 1.9 


0.22 ± 





02 


0.26 


0.26 


0.28 


0.24 


0.24 


0.24 


0.24 


K°/A 


1 


m+ > 1.9 


0.87 ± 





09 


1.33 


1.30 


1.38 


1.15 


1.16 


1.20 


1.20 


5/A 


1 


m± > 1.9 


0.17 ± 





01 


0.28 


0.27 


0.18 


0.17 


0.17 


0.18 


0.18 


S/A 


1 


m_i > 1.9 


0.38 ± 





04 


0.62 


0.64 


0.38 


0.38 


0.37 


0.30 


0.30 


n+n 
p/p 


1 


> 2.3 


1.7 ± 





9 


0.95* 


0.98* 


0.59* 


0.58* 


0.58 


0.52* 


0.52 


4 


Mid-rapidity 


11 ± 


2 




11.0 


11.2 


10.1 


10.6 


10.5 


7.96 


7.96 


A/p 


5 


4 TT 


1.2 ± 





3 


2.43 


2.50 


1.47 


1.44 


1.43 


1.15 


1.15 


h- 
p-p 


6 


4 TT 


4.3 ± 





3 


4.5 


4.4 


4.2 


4.1 


4.1 


3.6 


3.6 


h+-h- 
h+ + h- 


7 


4 TT 


0.124 ± 





014 


0.109 


0.114 


0.096 


0.103 


0.102 


0.092 


0.092 



^ S. Abatzis et al, WA85 Collaboration, Heavy Ion Physics 4, 79 (1996). 
^ S. Abatzis et al, WA85 Collaboration, Phys. Lett.B 347, 158 (1995). 
^ S. Abatzis et al, WA85 Collaboration, Phys. Lett.B 376, 251 (1996). 
* I.G.Bearden et al, NA44 Collaboration, Phys. Rev. C 57, 837 (1998). 
^ D.Rohrich for the NA35 Collaboration, Heavy Ion Physics 4, 71 (1996). 

® S-Ag value adopted here: T. Alber et al, NA35 Collaboration, Eur. Phys. J. C 2, 643 (1998); lhep-ex/971100l|l . 

^ A. lyono et al, EMU05 Collaboration, Nucl Phys. A 544, 455c (1992) and Y. Takahashi et al, EMU05 Collaboration, 

private communication. 



we have not enforced the strangeness conservation in the fits D (in bold face) and D' shown in table |T|. 
Enforcing strangeness conservation leads to fits denoted (D^, D^) with x^/dof > 3, and thus not offering 
an acceptable interpretation of the data. There indeed is also not a well defined minimum in the 
associated distributions shown by dashed lines in Fig.|l[ 

Thus the following simple reaction scenario emerges: in the particle evaporation-dissociation process 
at Tf = 143 lb 3 MeV there is predominant emission of s-carrying hadrons (see below) leading to a s-carrier 
rich residue that populate lower p± kaon and hyperon yields that remained not fully measured. It should 
be here stressed that the data we fit are very sparse in the small p± < 1 GeV domain, thus an imbalance 
between s and s we show in table |3|, simply should be read to mean that the opposite imbalance is to be 
expected for soft hadrons. A soft hadron excess arising from disintegration of s-quark enriched residue 
becomes particularly visible in the J7-abundance, since the abundance of these particles is small, while 
their production by strangeness rich source is probably enhanced, speaking here in relative terms. We 
note that such strangeness distillation has been foreseen to occur [17|, and specifically in this manner 
if QGP phase hadronizes at temperatures found in our fits of the data [^]. This dynamical distillation 
process stimulated several searches for strangeletts, i.e., (quasi) stable highly strange hadrons with more 
than three strange quarks [|^, . 
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Figure 1: x^/dof profile as function of temperature for the chemical non-equilibrium fits: without strang- 
eness conservation constraint D (thick solid line), D' (thin solid line), and with strangeness conservation 
constraint Dg (thick dashed line), Dg (thin dashed line). Thick lines exclude the three J7-data points. 



4 Properties of hadronic particle phase space 

Given the precise statistical information about the hadron phase space provided by the fit D, we can 
determine the specific content in energy, entropy, strangeness carried out by hadronic particles, see table |. 
For each fit we show the chemical freeze-out temperature Tf, the hadron phase space properties: specific 
(per baryon) energy, entropy, (anti)strangeness, strangeness imbalance, and the freeze-out pressure Pf 
and volume Vf. The specific s content is determined to be 0.90 it 0.04. We find that there is an excess 
of s-carrying hadrons carried away at higher -p^. The specific entropy content 48.2 it 3 of the best fit D, 
agrees well with the evaluation made earlier |jl^. This is so because the fit of predominantly high 
strange particle data is fully consistent with the 47r total multiplicity results. The pressure of the hadronic 
phase space Pi = 82±6MeV/fm'' has the magnitude of the vacuum confinement pressure. If the chemical 
freeze-out were the result of a disintegration after QGP expansion to the point of instability into volume 
fragmentation, this would be according to the Gibbs criterion the transformation pressure. In that case 
our fit amounts to the measurement of the confinement pressure (bag constant B\ The volume Vi we 
quote arises by assuming that the hadronic phase space comprises for the S-W collisions baryon number 
B = 120. It implies in the case of fit D a hadron phase space energy density ef = 0.43 it 0.4GeV/fm^. 
Both, Vf and ef have physical meaning only if the chemical freeze-out occurs nearly simultaneously within 
the entire body of the source. 

The energy per baryon found in emitted hadronic particles is remarkably consistent with the expec- 
tation based on kinematic energy content, should the deposition of energy and baryon number in the 



fireball be identical [21|. For the favored fit D, we find Ef/B = 9.05 it 0.5 GeV, to be compared with 
kinematic value 8.8 GeV obtained for central collision of Sulphur with a tube of matter in the larger 
target. Moreover, the different fits cross the kinematic value (dotted horizontal line in Fig. ^, bottom 
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Table 3: Tf and physical properties of the full hadron phase space characterized by the statistical 
parameters given in table ||. 



Fits 


Tf [MeV] 


Et/B 


Si/B 




(sf - Sf)/B 


Pf [GeV/fm^] 


Vi [fm^] 


A 


145 ± 3 


9.01 ± 0.50 


50.1 ± 3 


1.64 ± 0.06 


0.37 ± 0.02 


0.056 ± 0.005 


3352 ± 350 


B 


144 ± 2 


8.89 ± 0.50 


50.0 ± 3 


1.66 ± 0.06 


0.44 ± 0.02 


0.056 ± 0.005 


3343 ± 350 


C 


147 ± 2 


9.25 ± 0.50 


48.5 ± 3 


1.05 ± 0.05 


0.23 ± 0.02 


0.059 ± 0.005 


3529 ± 350 


D 


143 ± 3 


9.05 ± 0.50 


48.2 ± 3 


0.91 ± 0.04 


0.20 ± 0.02 


0.082 ± 0.006 


2524 ± 250 


D' 


144 ± 2 


9.07 ± 0.50 


48.3 ± 3 


0.91 ± 0.05 


0.20 ± 0.02 


0.082 ± 0.006 


2521 ± 250 


D, 


153 ± 3 


8.89 ± 0.50 


45.1 ± 3 


0.76 ± 0.04 


0* 


0.133 ± 0.008 


1520 ± 150 




153 ± 3 


8.87 ± 0.50 


45.1 ± 3 


0.76 ± 0.05 


0* 


0.134 ± 0.008 


1515 ± 150 



portion) very near to best fit temperatures. The computed values of Ef/B are shown in Fig.|2| for a given 
Tf for fit D (solid line, no O or strangeness conservation) and fit Dg (dashed line, no $7, with strangeness 
conservation), with all the other parameters obtained from a least square fit at given temperature. The 
locations of the best x^/dof are indicated by vertical lines. In second section from the bottom in Fig.|2|, 
we draw attention to another remarkable physical result of the fit D: the value As = 1.00ib0.02 suggests as 
the source of these particles a state symmetric between s and s quarks, naively expected in the deconfined 
phase. It is not an accident, and we could have easily found another result, since Ag varies as function of 
fitted Tf. 

The non-equilibrium parameters, shown in top two sections of Fig. ^, are opening the minimum in 
the x^/dof distribution shown in Fig.||, without a major dislocation in the parameter space as is seen in 
table |l[ The result of interest here is that both consistently differ from their equilibrium values 7s = 1, 
7q = 1 . The values 7s = 0.73 it 0.02 and 7q = 1.22 it 0.06 are compatible with the physics based 
expectations: the computed rate of strangeness production in the deconfined phase is suggestive of this 
result which is noticeably below the full phase space equilibration |10, 13 1. Fragmentation of remaining 
gluons at Tf = 143 MeV naturally leads to the value of 7q we obtain. As shown in the top half of 
Fig.^ for the entire range of freeze-out temperature of physical interest, the values of 7s and 7q vary 
somewhat, but consistently remain in the physical domain described above, and are similar for fits with 
and without strangeness conservation. This systematic clearly assures that these parameters are here 
to stay. Even a short view of Fig.^ shows that strangeness conservation condition, which introduces 
the greatest systematic uncertainty into this study, is capable to influence the value of As, but not the 
chemical non-equilibrium parameters shown in table ||, or other interesting properties of the hadron 
source displayed in table ^. 

The sensitivity of the computed specific negative hadron yield to the freeze-out temperature Tf is 
shown in Fig.|3[ For a given value of Tf all the other parameters follow values determined by least square 
fit. The solid line is for fit D, dashed line for strangeness conserving fit Dg. Vertical lines are indicating the 
location of the best fits, horizontal dotted lines show the range of the NA35-experimental result for S-Ag 
collisions. The result shown in Fig.^ makes it clear that it is impossible to consider chemical freeze-out 
to occur outside of the temperature window 135 < Tf < 155 MeV. Our use of all accessible experimental 
particle production data in S-Au/W/Pb reactions eliminates thus from further consideration the recent 
suggestion of hadronic particle production universality between e'^ + e~ , p + p, p + p, A + A reactions 
with chemical freeze-out at Tf ~ 180 MeV [^]. In our view, it is more natural to see earlier (i.e., higher 
Tf) chemical freeze-out in light collision systems, as compared to heavy ion collisions. We note that in |Q, 
as here, the chemical freeze-out conditions are studied, and that subsequent to establishment of particle 
abundances there could still be particle number conserving (elastic) evolution of the hadronic matter to 
the point of full decoupling at a lower temperature. During this evolution the final shape of the particle 
spectra is determined. Thus data fits based on details of particle spectra, in particular including the low 
m_|_, are geared to yield a thermal freeze-out temperature Tfh that is yet lower than our Tf [22|. 
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5 Impact of collective flow on freeze-out properties 

So far we have considered only the particle yields, and have not addressed particle spectra which, as 
mentioned in beginning of section requires that we allow for the Doppler-like blue shift of the particles 
emitted from a moving source at hadronization. To do this we need to explicitly add to the least square 
fit another parameter, the collective velocity Vc- While the integral over the entire phase space of the 
flow spectrum yields as many particles with and without flow, when acceptance cuts are present particles 
of different mass experience differing flow effects. We consider here the radial flow model, perhaps of the 
simplest of the complex flow cases possible [^], but it suffices to fully assess the impact of flow on our 
analysis. 

For given pair of values Tf and Vc, the resulting m± particle spectrum is obtained and analyzed using 
the shape and procedure employed by the experimental groups, and the inverse slope 'temperature' 
is determined for each particle j — since we have the experimental results for , we can indeed include 
these in our data fit. Vc is thus just one additional parameter and there are in principle several new data 
points available. However, in consideration of their similarity |^], with values overlapping within error, we 
decided to include in the chemical fit only one value Tg = 235ibl0, chosen near to the most precise Lambda 
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Figure 3: Variation of h~ /{p — p) as function of Tf for fits D (solid line) and strangeness conserving Dg 
(dashed line). Vertical lines point to best fit location, horizontal dotted lines brace the experimental 
result (see table |2|) . 

spectra. Once the chemical fit yields Tf and Vc, we check how other particles have fared. The resulting 
are in remarkable agreement with experiment, well beyond what was fitted: we find for kaons, lambdas 
and cascades the values = 215, 236 and 246 MeV respectively, which both in trend and value agrees 
with the K°, A, A, E and E WA85 results [||]: Tf " = 219±5, T^^ = 233±3, t]" = 232±7 , Tf = 244±12 
and Tf = 238 ± 16. 

Since the flow effect shifts particles of different mass differently into different domains of m± , y , it is 
not surprising that the inclusion of flow only impacts the phase space abundance parameters, beyond the 
errors of the fits. The fit D remains the best, it has now x^/dof=0.73, and the parameter values we gave 
in abstract. Fit D' has x^/dof=0.83 . In order to facilitate comparison with other work we note that to 
greater precision our fit D with flow yields T = 142.7 ± 2.1 MeV, Hb = 176 ± 3 and Vc = 0.486 ± 0.010 c. 
The situation with strangeness non-conservation remains unchanged, both fits Dg and D<j' have negligible 
confidence level with x^/dof~ 3.2. The strangeness imbalance in fits D, D' is the same as we have 
obtained without flow. Thus this imbalance is not result of flow effects, but is an intrinsic property 
of the evolution of the fireball. An interesting feature of with flow fit is that there is little correlation 
between now 6 parameters in the fit, in other words the flow velocity is a truly new degree of freedom in 
description of the experimental data. We checked that nearly the same flow velocity is found when we 
disregard in the fit the experimental inverse slope. It is for this reason that we present here the simple 
radial flow model, as within this scheme the inverse transverse slope of hadrons is correctly 'predicted' 
by the chemical freeze-out analysis with flow. We note that in just one aspect the fits with flow offer a 
new insight: the value of 7^ we obtained is compatible with unity. This value was noted already in the 
fit of the S-S collision results |^], and so far eluded the analysis of S-W/Au/Pb collisions. 

With 7g ~ 2 the analysis of excess entropy production is fully confirmed, and indeed in fit D 
with flow the entropy per baryon is 46.3 , within the bounds shown in table |3|. The energy content in 
the phase space is slightly different: the intrinsic energy per baryon at freeze-out is E/B = 8.55 GeV, 
which when multiplied with the factor 7^ reaches E/B = 9.78 GeV. The freeze-out pressure Pf = 0.123 
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GeV/fm^, roughly growing with a power of jq, has accordingly increased. 

As this discussion has shown, the analysis of particle abundances and spectra has to be combined in 
order to be certain that all the physical parameters acquire their true value. Even though there is still 
considerable uncertainty about other freeze-out flow effects, such as longitudinal flow (memory of the 
collision axis), the level of consistency and quality of our fit suggests that for the observables considered 
here these effects do not matter. Considering the quality of the data description obtained using the 
chemical non-equilibrium parameters, it is impossible to classify results presented here as accidental and 
likely to see further major revision. 

In conclusion, we have here shown that the hard > 1 GeV) strange particle production data, 
combined with the global hadron multiplicity (entropy), can be consistently interpreted within a picture 
of a hot hadronizing blob of matter governed by statistical parameters acquiring values expected if the 
source structure is that of deconfined QGP. We have shown that radial flow allows to account exactly for 
the difference between freeze-out temperature and the observed spectral shape, and allows full description 
of inverse slope of m± strange baryon and antibaryon spectra. We find a highly significant fit to data with 
a source characterized by As = 1, specific baryon energy E'/i? ~ 9 GeV, high specific entropy S/B ~ 50, 
and chemical freeze-out temperature Tf = 143 MeV . The dense blob of matter was expanding with 
surface velocity Vc — 0.49 c. The near equilibrium abundance of strange quarks (7s — 1, including flow), 
and the over-abundance of light quarks (7q ~ 2), is pointing to a deconflned, fragmenting quark-gluon 
fireball as the direct particle emission source. At freeze-out we find the particle phase space pressure 
Pf ~ O.lGeV/fm^, which maybe interpreted as the pressure needed to overcome the external vacuum 
pressure (bag constant B). Strangeness non-conservation at high p± seen in our fit favors an evolution 
scenario in which the process of hadronization into hard particles enriches by distillation the fireball 
remnant with strangeness content. Given an excess of s emission at high pj_ of 0.2 strange particles per 
baryon, we expect at small p± a non-negligible 20% excess of s-carriers over s-carriers, potentially helping 
the formation of strange quark matter nuggets. The lessons learned in this study will likely prove helpful 
in the forthcoming interpretation of the Pb-Pb data, and the future RHIC-BNL data analysis. 
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